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ABSTRACT 

Recent submillimetre surveys have revealed a population of dusty, high redshift sources 
of great cosmological significance for understanding dust-enshrouded star formation 
in distant galaxies, and for determining the origin of the far-IR background. In this 
paper, we analyze nine rich cluster fields mapped at 850 /im and 450 /im with the 
SCUBA array on the James Clerk Maxwell telescope. Lensing models of the clusters 
are developed in order to derive accurate source counts for our sample. VLA maps of 
the same clusters are used to help constrain the redshift distribution of our SCUBA 
detections. Implications for high redshift galaxies and for the far-IR background are 
discussed. We also provide limits on distributed dust produced by cooling flows in 
these clusters. 
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1 INTRODUCTION 

The submillimetre waveband has recently become an invalu- 
. 1^ ' able tool for investigating the properties of high redshift 
^^ , galaxies. An early result using the Sub-millimetre Common 
'j_| ■ User Bolometer Array (SCUBA - Holland et al. 1999) by 
Cu ^ Small, Ivison and Blain (1997) showed that a much larger 
population of dusty, high redshift galaxies existed than 
previously thought. The existence of these dusty, sub-mm 
galaxies at relatively high redshifts have now been confirmed 
in several surveys (Hughes et al. 1998; Barger et al. 1998, 
1999a; Bales et al. 1999; Lilly et al. 1999a; Blain et al. 1999). 
Such sub-mm sources could be of great cosmological signif- 
icance - they may comprise a substantial fraction of the 
Far-Infrared Background (FIRB - Puget et al. 1996, Hauser 
et al. 1998, Lagache et al. 1999). Additionally, they could 
be an important source of CMB anisotropy at arcsecond 
scales (Scott & White 1999), and certainly the dominant 
source of sky fluctuations at ~ 1mm, at least out of the 
galactic plane (Hughes et al. 1998, Borys, Chapman & Scott 
1999). Our understanding of the star- formation history of 
galaxies depends critically on the properties of these sub- 
mm sources: this population is responsible for a substantial 
fraction of the total star formation output in the early Uni- 
verse and holds important information about how galaxies 
formed and evolved (e.g. Blain & Longair 1993, Guiderdoni 
et al. 1997, Blain et al. 1998). They appear to be the high 
redshift counterparts of well-studied nearby ultra-luminous 



IR galaxies (Sanders & Mirabel 1996), and can be plausi- 
bly intrepretted as the formation sites of galactic spheroids 
through mergers of gaseous disk systems (e.g. Sanders 1999, 
Lilly et al. 1999a). 

Much work has gone into identifying just what types of 
galaxies make up this sub-mm bright population, and what 
mechanisms may be responsible for their rest-frame far-IR 
emission (Smail et al. 1998; Hughes et al. 1998; Ivison et 
al. 1998b; Lilly et al. 1999b; Frayer et al. 1998; Frayer et 
al. 1999; Frayer et al. 2000; Chapman et al. 1999; Ivison et 
al. 2000; Chapman et al. 2000; Bertoldi et al. 2000; Gear 
et al. 2000). However, the numbers of sub-mm source de- 
tections are still relatively small. In an effort to increase the 
statistics and broaden the source count baseline, we have un- 
dertaken a sub-mm survey including both blank sky (Borys 
et al. in preparation) and cluster fields. Our cluster sample 
is complementary to that of Smail, Ivison, and Blain (1997, 
1998) who pioneered this approach using SCUBA on the 
JCMT. These cluster fields benefit from the lensing by the 
foreground cluster mass at z<l, enhancing sensitivity to 
the dusty galaxies at higher redshifts. 

Our cluster sample consists of nine clusters at redshifts 
0.2-0.8, originally selected to search for both distributed 
cool dust (cooling fiows - see Edge et al. 1999) and the 
Sunyaev-Zel'dovich (SZ) effect (Birkinshaw 1999; Andreani 
et al. 1999). Although we found no convincing evidence for 
distributed emission, our SCUBA maps have revealed sig- 
nificant source detections, as expected from the counts from 
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other surveys (e.g. Blain et al. 1999, Barger et al. 1999b). 
This further demonstrates that sub-mm surveys provide an 
efficient means of identifying a population of dusty AGN 
and/or star-forming galaxies, possibly at very high red- 
shifts. In addition, a few weak sources appear coincident 
with radio-bright central galaxies in the clusters, possibly 
associated with the strong cooling flows suggested by the x- 
ray profiles. In this paper, we describe possible sources found 
in the cluster fields, and present the resulting source counts. 
Existing VLA maps are used to place limits on the redshift 
distribution of the sample. We conclude by discussing the 
implications for high redshift galaxies and the FIRE. 



same cluster field, and sky rotation, ensure that (in most 
cases) each point of the map is covered by many bolometers. 
Even with the inverse variance weighting, excessively bad 
pixels still appeared to degrade the map, and were therefore 
clipped from the data prior to rebinning into the final maps. 
The data were then calibrated against standard planetary 
or compact Hll region sources, observed during the same 
nights as the clusters. Uncertainties at 850 /im and 1350 /im 
are estimated to be around 10 per cent, while at 450 pim 
they could be as high as 25 per cent. Note that calibration 
uncertainties do not affect the S/N level for our detections, 
only the overall flux densities reported. 



2 OBSERVATIONS 

2.1 JCMT data 

Nine clusters were observed with the SCUBA instrument on 
the James Clerk Maxwell Telescope, resulting in 17 signifi- 
cant source detections greater than 3cr. During 4 observing 
runs throughout 1998 and 1999, we operated the 91 element 
Short-wave array at 450 /im, and the 37 element Long- wave 
array at 850 /im simultaneously, in JIGGLE mode, and ad- 
ditionally the single photometry pixel at 1350 /im for some 
sources, giving half-power beam widths of 7.5, 14.7, and 21.0 
arcsec respectively. Additionally, 850/450 ^m PHOTOME- 
TRY mode foUowup on 3 of the sources was performed, con- 
firming their existence and flux densities. The central pixel 
of SCUBA was fixed on the source, as defined by the centroid 
of the ROSAT imaging (Crawford et al. 1999). For mapping, 
the standard 64-point jiggle pattern was employed to fully 
sample the 450 and 850 /im arrays. The effective integra- 
tion times (time spent on and off source derived from the 
number of 64 point integrations obtained, which excludes 
overheads) are presented in Table 1 for each cluster. For 
1350 /xm photometry we used the 9-point jiggle pattern to 
reduce the impact of pointing errors by averaging the source 
signal over a slightly larger area than the beam, resulting in 
greater photometric accuracy. Whilst jiggling, the secondary 
was chopped at 7.8125 Hz with chop throws between 40 and 
120 arcsec in azimuth, with smaller throws used under worse 
atmospheric conditions or when the projected core radius for 
a cluster was smaller. 

Pointing was checked hourly on blazars and a sky-dip 
was performed to measure the atmospheric opacity directly. 
The RMS pointing errors were below 2 arcsec, while the 
average atmospheric zenith opacities at 450 pm , 850 fim and 
1350 pm were fairly stable and generally quite good with rsso 
ranging from 0.15 to 0.36. Some short time-scale variations, 
presumably due to the atmosphere, caused some parts of 
the data-set to be noisier, and these noisiest sections were 
excluded. 

The data were reduced using the Starlink package 
SURF (Scuba User Reduction Facility, Jenness et al. 1998). 
Spikes were first carefully rejected from the data, followed 
by correction for atmospheric opacity and sky subtraction 
using the median of all the array pixels, except for obviously 
bad pixels (the 1350 ^m pixel currently has no provision for 
subtracting sky variations using the other wavelength pix- 
els). We then weighted each pixel by its timestream inverse 
variance, relative to the central pixel. Multiple scans on the 



2.2 Archive VLA data 

VLA maps at 4.9 GHz in the C configuration were obtained 
for five clusters in our sample courtesy of A. Edge. The effec- 
tive resolution in the images is 5 arcsec. More details of the 
reduction and analysis of similar maps are given in Edge et 
al. (1999a). Typical la sensitivities lie around 60/iJy. In ad- 
dition, an archive 1.4 GHz map for Abell2219 was obtained 
courtesy of F. Owen, reaching a la limit of ~39/iJy per 2 
arcsec beam. 

For the remaining clusters, VLA data presented in 
Stocke et al. (1999) at 1.4 GHz in the B configuration, was 
used to set limits and identify counterpart positions. Typical 
1.4 GHz sensitivities in these maps are 65 ^Jy RMS. 



3 SOURCE CATALOGUE 

The 850 pim SCUBA maps for the clusters showing signifi- 
cant source detections are presented in Figs. 1 and 2. Any 
real objects in our maps are expected to be unresolved and 
will typically appear as a positive point spread function with 
FHWM --^ 15 arcsecs at 850 /im . For most of our maps, we 
do not expect to see negative signatures for our sources, as 
the individual scans coadded in the map have different chop 
throws and directions, diluting the effect of off-beams. We 
also chopped out of the field in most cases. The edges of 
the maps incur a rapid rise in noise level due to incomplete 
sampling of the sky for the outer bolometers of the SCUBA 
array. Eales et al. (1999) simulated the incomplete sampling 
at the field edge and found that any spurious sources in 
Monte Carlo simulations tended to appear in this region. 
We have therefore clipped the 14 arcsec edge regions in our 
maps where the noise level increases by a factor ^ 2 times 
that at the field center. 

As the signal-to-noise ratio (S/N) is low for individual 
pixels, we use a template fitting method for source detection 
to achieve a maximum gain in S/N, similar to that used in 
Eales et al. (1999). The cluster maps are convolved with the 
beam map obtained by imaging Mars, and all peaks greater 
than 3(7 are extracted as potential sources. We checked the 
convolved maps carefully against the raw data to ensure 
that no residual bad pixels were identified as sources. We 
also note that there are no strong negative sources in our 
maps. 

Sensitivity limits for deep SCUBA integrations have 
now been fairly well characterized (Small et al. 1997, Hughes 
et al. 1998, Eales et al. 1998, Barger et al. 1999b). But at 
these low S/N levels, we must be careful to consider possible 
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Figure 1. Cluster field images at 850 ^in, convolved with the SCUBA 14.7 arcsec beam. Dark shading is positive fiux density, 
while lighter shading is negative. All contour levels indicate the number of sigma above the background noise estimated from the 
central regions of the beam-convolved images, starting at 2cr, with intervals of la. These noise levels correspond approximately to 
the tabulated values (Tables 1&:2), which were derived from the raw data. Top row: CI 0016+16 (left); MS 0451-03 (right). Bottom 
row: Abell520 (left); Zwicky3146 (right). 



systematic noise effects. There are currently five identified 
components which contribute to the overall SCUBA noise 
level: (i) photon noise from the sky; (ii) telescope and instru- 
mental noise; (iii) blanked-off detector Johnson, phonon and 
amplifier noise; (iv) excess photon noise from the optics; (v) 
excess radiative load in the dewar. The noise terms add in 
quadrature, and the estimate for the relative contributions 
at the average sky emissivity for our survey ((rgso) ~0.2) 
are respectively 29, 25, 17, 13, and 16 per cent (W. Hol- 
land private communication). Thus sky noise is the largest 
source, but does not dominate the others. Nevertheless, we 
expect the noise to integrate down with time. This will not 
be true of one additional source of noise: 'confusion noise' 
arising from fluctuations due to undetected sources. How- 
ever, the level is essentially negligible for these observations 



- Blain et al. (1998) quote a value of aconf ~ 0.44 mjy, derived 
from their source counts, while we find that any reasonable 
extrapolation of the counts (together with the constraint of 
not over-producing the FIRE) leads to a variance of no more 
than 0.5 mJy. We can therefore safely ignore confusion noise. 
However, it is worth pointing out that lensing by the clus- 
ters effectively allows us to probe below the confusion limit 
of say 3(Tconf . 

With transient noise and bad pixels from scan to scan, 
our confidence in many of the detections is reinforced by 
dividing the data into two parts and seeing evidence of the 
sources (at the > 2a level) appear in both. Multi- wavelength 
detections also substantially strengthen our confidence in a 
source. Although the 450 /im channel is typically not sensi- 
tive enough to detect these objects, especially those lying at 



© 0000 RAS, MNRAS 000, 000-000 



4 S. C. Chapman et al. 



—\ — I — I — I — \ — I — r — r — T — \ — r-\ — } — rn — [— 




22°22'00" — 



r-| — 1 — 1 — \ — \ — I — I — r — \ — \ — I — 1 — T — I — r 



22°2r00" 



22°21J'00" 




22°19'00" 



1^ 



t32°09'00" 



t32°08'00" 



- t32°fir[ 



14^57"20= 18= 16= 



- ■ ' I 



12= 



I 




I , , , , I 



_L 



16"40"23= 26= 24= 22= 20= 18= 16= 14= 12= 



17^22"32= 30= 28= 26= 24= 22= 



Figure 2. SCUBA maps at 850 ^m continued. Contours arc the same as Fig. 1. Top row: MS 1054-03 (left); MS 1455+22 (right). 
Bottom row: Abcll2219 (left); AbeU2261 (right). 



higher redshifts, many more of them may be observable at 
1350 pim with deeper photometry than we have currently 
obtained. Three sources have been subsequently verified 
through PHOTOMETRY foUowup observations (MS 0451- 
A, Abell2261-A, and MS1455-A). These are discussed in 
the individual source section (4.2). 

Empirically, the localized noise levels are seen to inte- 
grate down in our maps roughly as i~^ " (see also Ivison et 
al. 1998a). We note that in a typical 850 ^m SCUBA field of 
5.2 square arcmin with ~ 100 independent beams we expect 
only ~0.1 spurious 3a peaks from Gaussian noise, and and 
<C 1 spurious 4a peak from our whole survey. 

The completeness of a set of detections has typically 
been established using synthetic source recovery techniques 
in Monte Carlo simulations (Small et al. 1997, Hughes et al. 
1998). Our 3a and source catalogues corresponds approxi- 
mately to the 80 per cent completeness level found in Small 



et al. (1997). We estimate that our 4a catalogue is more 
than 90 per cent complete. 

Our sub-mm observations of the clusters are presented 
in Tables 1 and 2. Table 1 lists the clusters, coordinates, and 
redshifts, along with total exposure time (based on the num- 
ber of integrations). The last four columns present the aver- 
age 850 //m and 450 /im sensitivities achieved after subtract- 
ing sources > 3a from the field, and the number of sources 
detected at greater than the 3a and 4a levels. Table 2 lists 
the properties of the detected sources, with coordinates and 
observed flux densities in each of the three sub-mm wave- 
length bands. These sub-mm detections and upper limits 
(upper limits are given as 95 per cent Bayesian confidence re- 
gions, assuming Gaussian errors and neglecting the negative 
fiux region) at 1350, 850, and 450 pim provide rough photo- 
metric estimates of the source redshifts (see e.g. Hughes et 
al. 1998). However, the addition of radio data significantly 
increases confidence in these photometric redshifts. The final 
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two columns list the VLA radio flux limits, and estimates of 
the redshift from the sub-mm and radio data, as described 
below. 



4 PROPERTIES OF INDIVIDUAL SOURCES 

4.1 VLA radio counterparts 

Radio counterparts to sub-mm sources are expected due to 
the synchrotron emission resulting from either an AGN or 
shocked gas in supernovae remnants from a starburst (e.g. 
Richards et al. 1999). Given the small number density of ra- 
dio emitting objects, the chances of finding a 1.4 GHz VLA 
object of > 200 /iJy in our error circle by chance is much less 
than 1 per cent (Richards et al. 1998). By contrast, identi- 
fying optical counterparts to the SCUBA-selected galaxies 
is a non-trivial matter for two main reasons: 1) the beam- 
size at 850 /im (the optimal wavelength for these studies) is 
~ 15 arcsec, with pointing errors for the telescope of order 
2arcsec; 2) the large, negative K-corrections (i.e. increase 
in flux density as the objects are redshifted) of dusty star- 
forming galaxies at these wavelengths imply that sub-mm 
observations can detect such objects at z > 1 in an almost 
distance-independent manner. Thus several candidate opti- 
cal galaxies are often present within the uncertainty of the 
sub-mm detection, along with the possibility that the actual 
counterpart is at much higher redshift and undetectable with 
current optical imaging. In this paper we restrict our atten- 
tion to radio counterparts, and defer detailed optical analysis 
of the sub-mm source positions to a separate paper. 

CariUi & Yun (1999, 2000 - CY) recently demonstrated 
that crude redshift estimates for distant dusty galaxies could 
be obtained in the absence of an optical counterpart using 
the spectral index between the sub-mm (850 /im) and radio 
(1.4 GHz) wavebands (af^^), with a mean galaxy model pa- 
rameterized as z = 0.050-0.308Q-f 12.40^ -23.00^ + 14.90". 
The model SEDs used by CY are based largely on lower 
redshift observations, however they are able to successfully 
describe the few high-redshift galaxies which have multi- 
wavelength radio and sub-mm observations available. Follow 
up studies by Blain (1999) have shown that if lower dust 
temperatures are adopted for the sub-mm population than 
are seen in the local sources used in CY's models, then the 
allowed redshifts are slightly lower for a given value of af^i" . 
The isolated cases which do not fall within the CY model 
predictions are known to have strong radio-loud AGN com- 
ponents. The important point is that the derived redshift 
provides a robust lower limit, irrespective of the nature of 
the emission mechanism, AGN or starburst, since an AGN 
contribution to the radio emission flattens the radio spec- 
trum. The modest scatter between the models in CY sug- 
gests that the af^f technique currently provides the most 
useful limit on the redshifts of sub-mm galaxies, in the ab- 
sence of an optical counterpart. 

Small et al. (2000) recently applied the CY analysis to 
deep radio observations of the complete sample of SCUBA- 
selected galaxies from their cluster lens survey (Small et 
al. 1998, Ivison et al. 2000), revealing a higher redshift dis- 
tribution than previously expected from optical follow-up 
work. Using existing VLA data for our SCUBA cluster sam- 
ple, we can thus constrain the redshift distribution of our 



population and set the stage for optical/ncar-IR follow-up 
identifications. 

We searched for radio counterparts around the nomi- 
nal positions of the sub-mm sources based on the SCUBA 
astrometry. Positional errors for SCUBA sources from the 
literature which were clearly identified in radio and optical 
wavelengths (e.g. Ivison et al. 1998b, Frayer et al. 1999) were 
~ 6 arcsec for 4cr detections and ~ 8 arcsec for 3a detections, 
and we use similar error circles for our analysis. Radio lim- 
its at 1.4 GHz for sources in A 520, Zw3146, MS 1455, and 
A 2261, have been extrapolated from the value at 4.9 GHz, 
using a spectral index a = 0.8 (i.e. S oc i^""®) appropriate for 
star formation driven synchrotron emission. This provides a 
reasonable estimate of the source redshift using the CY re- 
lation (since a brighter radio counterpart leads to a lower 
redshift estimate, while an AGN component would flatten 
the spectrum giving rise to a higher redshift). The 1.4 GHz 
VLA data from Stocke et al. (1999) are used to derive limits 
for CI 0016, MS 0451, and MS 1054, and an archive 1.4 GHz 
map of A 2219 (courtesy of F. Owen) was used to measure 
fluxes directly. In 5 cases, there were radio sources > 3a 
lying within the error circle of the sub-mm beam. One of 
these may be a central cluster galaxy and is also discussed 
in the following section. The 15 arcsec by 15 arcsec region 
surrounding the sub-mm sources are shown in Fig. 3, with 
the apparent radio flux densities or limits listed in Table 2. 

We use the CY relation (which is consistent with dust 
temperature Td — 50 K) to generate these lower limits to the 
source redshifts. We note that the actual dust temperature 
in the sub-mm sources may be lower, bringing down the 
redshifts, or of course higher, thereby raising the inferred 
redshift. In general the constrained quantity is Td/(1 + z). 



4.2 Source characteristics and evidence for cluster 
members 

We now discuss each cluster field in detail: 
C10016-hl6 

The data for this cluster show a general positive excess at 
the ~ 2a level, with 2 peaks barely reaching the 3a level in 
the convolved map (sources A and B). There are no 1.4 GHz 
VLA counterparts down to the 200 /iJy level, although the 
bright 2.7 mjy radio souce from Stocke et al. (1999) lies 
on a clear sub-mm peak which reaches 2.6a significance in 
the convolved map. The generally positive excess could be 
caused by a Sunyaev-Zel'dovich increment (core radius ~ 2 
arcmin - Carlstrom et al. 1999), which is too noisy in this 
map to detect with any confidence. However it may have 
some effect on the detectability of other sources in this field. 
MS 0451-03 

An unresolved source (B: SMMJ 04542-0301) is found to 
the West in this cluster, and there is north-south extended 
positive sub-mm emission over the x-ray centroid region, 
reaching 4a in the convolved map (A: SMMJ 04541-0302). 
An additional 850 /im photometry observation centered on 
the elongated source reveals a detection of 10.1±3.0mJy, 
which is consistent with the peak flux densities found along 
the source as well as the positional accuracy of the SCUBA 
map. This central source is offset by 20 arcsecs to the east 
from the brightest cluster galaxy (BCG), and 20 arcsecs to 
the west from a 1.8 mJy VLA source (Stocke et al. 1999). 
This detection could possibly represent cool dust associated 
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(J2000) 



Tabic 1 
Sub-mm observations for the clusters in our sample 

eb 



Doc." 
{,72000) 



Exposure time 
(ks) 



(mjy) 



(mJy) 3cr 4cr 



Cl 0016+16 


00 18 33.2 


+ 16 26 18 


0.541 


20.5 


3.1 


36 


2 





MS 0451-03 


04 54 13.4 


-03 01 47 


0.550 


14.1 


4.2 


181 


2 


2 


Abell520 


04 54 19.0 


+02 56 49 


0.203 


7.7 


5.0 


66 


2 


1 


Zwicky 3146 


10 23 39.6 


+04 11 10 


0.291 


12.8 


2.6 


25 


4 


1 


MS 1054-03 


10 57 00.2 


-03 37 27 


0.833 


7.7 


4.3 


103 


1 





MS 1455-1-22 


14 57 15.1 


+22 20 35 


0.259 


15.4 


2.1 


22 


1 


1 


Abell2163 


16 15 34.1 


-06 07 26 


0.201 


7.7 


5.8 


44 








Abell2219 


16 40 20.5 


+46 42 59 


0.228 


15.4 


1.7 


8 


4 


1 


Abell2261 


17 22 24.1 


+32 07 45 


0.224 


7.7 


3.8 


105 


1 


1 


Total 














17 


7 



^ x-ray centre — ROSAT archive. 

Measured RMS over the central 60 arcsec of the field after subtracting a scaled SCUBA beam profile from any identified sources in the region. 



Tabic 2 
Details of individual sources 



Source 


qb 


qb 
■^850Mm 


qb 


^l.i GHz 


4.. 






(mJy) 


(mJy) 


(mJy) 


(pJy) 


C10016 












A SMMJ 00186+ 1626 




9.9+3.3 


7 + 34 


<200 


>1.8 


B SMMJ 00186+1625 




9.8+3.2 


- 4 + 36 


<200 


>1.8 


MS 0451 












A SMMJ 04541-0302 




16.8+4.2 


24 + 181 


<200 


>2.2 


B SMMJ 04542-03018 




19.1+4.2 


-61 + 177 


<200 


>2.3 


Abell520: 












A SMMJ 04543+0257 




33.0+5.0 


-39 + 67 


<420 


>2.2 


B SMMJ 04543+0256 




16.2+5.1 


13 + 70 


<410 


>1.6 


Zwicky 3146 












A SMMJ 10237+0411 


1.7 + 3.9 


9.3+2.5 


15 + 23 


482/2722° 


1.1(0.3)° 


B SMMJ 10237+0410 




9.2+2.5 


-2 + 25 


543 


1.1 


C SMMJ 10237+0412 




9.4+2.6 


6 + 27 


<450 


>1.2 


D SMMJ 10236+0412 




11.2+2.6 


23 + 27 


467 


1.3 


MS 1054 












A SMMJ 10571-0337 


15.2+5.4 


14.8+4.3 


12 + 103 


<200/1600 


2.1(0.8) 


MS 1455: 












A SMMJ14573+2220S 




10.1+2.2 


5 + 22 


<294 


>1.5 


Abell2219 












A SMMJ 16403+46440 


4.3+2.0 


10.6+1.7 


34 + 11 


<220 


>1.8 


B SMMJ 16403+46437 




5.8+1.7 


- 2 + 8 


593 


0.8 


C SMMJ 16404+4643 




<5.7 


26 + 8 


280/9100 


0.9(0.2) 


D SMMJ 16404+4644 




6.3+2.0 


n/a* 


1270 


0.4 


Abell2261 












A SMMJ 17223+3207S 




17.6+3.9 


-56 + 105 


<452 


>1.6 



" Positional errors for SCUBA sources with identifications at radio and optical wavelengths (e.g. Ivison et al. 1998b, Fraycr et al. 1999) are ^ 6 arcsec for 
4cr detections and ^ 8 arcsec for 3cr detections. 

Estimated flux densities with Icr errors. All non-detection upper limits are given as 95 per cent Bayesian confldence regions, assuming Gaussian errors and 
neglecting the negative flux region 

° Bayesian 95 per cent limits at 1,4 GHz for sources in A 520, Zvi' 3146, MS 1455, A 2219, and A 2261, have been extrapolated from the value at 4.9 GHz, using 
a spectral index CK — 0.8 appropriate for star formation driven synchrotron emission. For the data from Stocke et al. (1999), only 3(T upper limits are given. 

Redshift estimate using Carilli & Yun (1999) radio/far-IR correlation. 
® In eases ^vhere there is ambiguity bet^veen the central cluster galaxy and a background galaxy, the radio flux of the cluster galaxy is also listed, and the 
cluster redshift is given in brackets. 

This source falls at the very edge of the short wavelength array. 
s These sources vi'cre verified with 850/450 ^m photometry. See text for details. 



with a cooling flow or could be associated with the Sunyaev- 
Zel'dovich effect - in either case the emission would have to 
be very strong an isolated to a small region of the cluster. 
For this particular cluster, the in-field chop reveals signs of 
the negative beams. However, another map of this cluster, 
observed with a large chop for Sunyaev-Zel'dovich increment 
measurement (Halpern et al. in preparation) suggests that 
additional source structure in the central region would sig- 
nificantly distort a clear offbeam signature. The extended 
nature of this central source (also verified in our indepen- 
dent map) is most likely to be a blend of two or more sources. 



For the purposes of Table 2 we treat this as a single source. 
Positions and fiux density estimates for such blends should 
be treated with some caution until higher resolution images 
are available. 
AbeI1520 

A large gradient appears to be present across the field, as 
evidenced by the increased positive excursions in the lower 
halves of both the 450 and 850 /im images and negative 
excursions in the upper half. This could either be real or 
due to some residual atmospheric signal - we could not 
find a convincing way to tell from our data, although we 
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Figure 3. VLA counterparts to our sub-mm detected galaxies. We siiow VLA sources for cases wliere a > 3cr VLA detection is 
seen relatively nearby the SCUBA position. The radio contours are at 3, 5, 10, 20, 50, 100 and 200(T in all cases. The circles are 
drawn at the positions of sources (see Table 2) from our > Ser SCUBA source catalogue, with 10 arcsec diameter as an estimate 
of the positional uncertainty. The radio maps are at 4.9 GHz except for Abcll2219-C which is at 1.4 GHz. The other 12 SCUBA 
sources have only VLA upper limits, as listed in Table 2. 



suspect the latter explanation. An apparent 4a source (B: 
SMMJ 04543+0256) lying to the south of the primary source 
(A: SMMJ 04543+0257), is ~3a above the noise even once 
a linear gradient is fit to the plane and subtracted. The 
flux density presented in Table 2 for this source reflects the 
subtraction of this linear gradient component. There is no 
detection of either the BCG or the bright radio source lying 
to the east in the archive VLA map. 
Zwicky 3146 

The central source in Zw3146 (A; SMMJ 10237+0411) ap- 
pears to be extended east-west, and may be another blend. 
The VLA map reveals that the western extent of the SCUBA 
source is associated with the BCG, with fairly bright radio 
emission (5'i.4ghz =0.92mJy). However, the eastern part of 
the SCUBA source shows no radio emission (< 210 /iJy), and 
the CY analysis suggests a high redshift background source. 
Sources B and C in this cluster field also have weak 4.9 GHz 
detections in the archive VLA map. 
MS 1054-03 

The source in the MS 1054 field (A: SMMJ 10571-0337), 
has a 1350/850 jj,m fiux density ratio near unity. As ther- 
mal dust emission alone would produce a steep oc ly^'' spec- 
trum, this may indicate that the object has a flat-spectrum 
AGN radio contribution at 1350 /im, and perhaps at 850 /im 
as well. A bright VLA counterpart, offset by 15 arcsec 
from the sub-mm peak, has been identified by Condon et 
al. (1998) as a cluster AGN radio source, and the CY re- 
lationship is consistent within errors with the source lying 



in the cluster at z~0.8. However, such a large radio/sub- 
mm offset is unlikely (see for example Richards et al. 1999), 
and SMMJ 10571—0337 could easily be a background galaxy 
lensed by the cluster. Within errors, the 1350/850 /im ratio 
is not inconsistent with a higher redshift source. This source 
has been verified in another SCUBA map of this cluster, al- 
though at a slightly lower flux level (P. van der Werf, private 
communication) . 
MS 1455+22 

There is some evidence for 850 /im SCUBA flux associated 
with the BCG in the raw MS 1455 maps. However, a some- 
what noisy central pixel during the observations of this clus- 
ter skews the peak proflles sufficiently that our source find- 
ing routine does not identify this as a detection (the source 
appears at only 2.6a in the final convolved map). We do not 
include this as a source in our catalogue until it can be veri- 
fied. However, the proximity of the sub-mm peak to 2 bright 
radio sources (Fig. 4), and a possible CO (1-0) detection for 
the BCG (Edge et al. in preparation) lends credence to this 
putative source. The bright sub-mm source to the south (A: 
SMMJ 14573+2220) has no VLA counterpart in the archive 
maps. An additional 850 /xm photometry observation cen- 
tered on source-A reveals a detection of 6.8+1.8 mjy, which 
is consistent within the relative errors between the mapping 
and photometry measurements, as well as the extended na- 
ture of this source and positional uncertainties. 
AbeI12163 
No sources were found in this cluster field down to a 3a 
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sensitivity of ITmJy (corresponding roughly to 6niJy after 
correction for lensing). This cluster had the poorest sensi- 
tivity, and we include it here for completeness only. 
AbeI12219 

This map has by far the lowest 450 /im rms, and as a re- 
sult there are, in addition to the 850 /xm detections, two 
450 pm detections at the 3a level. The 1350/850 fim and 
450/850 pm ratios and lack of a 4.9 GHz detection for source 
A (SMMJ 16403+46440) are consistent with the high red- 
shift predicted from CY. This is the only source in A 2219 
that does not have at least positive flux density in the VLA 
archive map. The central source (C: SMMJ 16404-1-4643), is 
~ 2.2(7 at 850 /im but has a fairly clear detection at 450 /im. 
It does not quite align with the central galaxy radio emission 
(see Fig. 4), but within errors, it is possible that the source is 
a cluster member. Indeed the large implied 450/850 /im flux 
density ratio is indicative of a rather low redshift source. 
In the archive 1.4 GHz map, however, there is a 5a source 
lying near the SCUBA centroid, offset ~ 7 arcsec from the 
BCG. Source D (SMMJ 16404-^4644) lies near the edge of 
the frame (and off the 450 jj,m frame entirely), and is only 3a 
above the local noise level. However the centroid is within 

5 arcsec of an ISO 15 /im source (Barvainis, Antonucci & 
Helou 1999), with which the sub-mm source may be associ- 
ated. 

Abell2261 

Although the bright (> 4a) source to the west (A: 
SMMJ 17223+3207) does not correspond with any radio 
emission, the central 2.5a peak in the SCUBA map is iden- 
tified with a radio source coincident with the optical BCG 
(Fig. 4). An additional 850 /im photometry observation cen- 
tered on source-A reveals a detection of 13.6+2.8 mjy, which 
is consistent within the relative errors between the mapping 
and photometry measurements, as well as the extended na- 
ture of this source and positional uncertainties. 

4.3 Cooling flows and centrally concentrated dust 
emission 

The high central densities of the x-ray gas in many lumi- 
nous clusters lead to predicted cooling times shorter than 
the Hubble time (e.g. Fabian 1994). This cooling gas, which 
can not support the pressure exerted by the overlying hot 
gas, is expected to flow into the cluster centre. Supporting 
evidence for such cooling flows is now available in the form 
of inverted x-ray temperature profiles (e.g. Buote 1999), and 
x-ray emission lines show gas that has lost most of its ther- 
mal energy. The fate of this cool gas remains controversial 
(e.g. Braine & Dupraz 1994) as no sink has yet been clearly 
identified. 

Observations of the central cluster galaxies in cooling- 
flow clusters indicate that the cooling flows are not form- 
ing large numbers of stars. On the other hand, such stars 
could be hidden if the initial mass function is biased towards 
low masses, making the integrated stellar spectrum difficult 
to detect (e.g. Fabian, Nulsen & Canizares 1982; Mathews 

6 Brighenti 1999). However, the emission-line nebulosity 
around the central galaxies in many cooling-flow clusters, 
seen on scales of up to 100 kpc, have emission-line ratios 
suggestive of large amounts of dust in the central regions 
- this dust could obscure any stars formed there (Hansen, 
Jorgensen & Norgaard-Nielsen 1995; Allen et al. 1995). HST 



imaging of central cluster galaxies (e.g. for Abell 1795 - Mc- 
Namara et al. 1996; Pinkney et al. 1996, for AbeU2597 - 
Koekemoer et al. 1999) has also directly detected dust lanes 
in some cases. At present the total amount or temperature 
of this dust is difficult to estimate from the scant observa- 
tional evidence. Whether the observed dust originate from 
the cooled gas clouds themselves or from the on-going star 
formation remains unknown. 

SCUBA data may help here by probing emission from 
cool dust in or around the central galaxies in clusters. The 
implications of detected cool dust in the BCG of massive 
clusters for the interpretation of star formation in cooling- 
flow galaxies, has been discussed in Annis & Jewitt (1993) 
and Edge et al. (1999) Although there are no highly sig- 
nificant detections of the BCGs in our sample, 3 of the 
clusters show a marginal positive flux density at 850 /xm 
(~2.5a) near the BCG. In Abell 2219, a 3.5a detection at 
450 ^im (2.3a at 850 /im) lies very near a radio luminous 
BCG (Fig. 4). In addition, the cluster MS 0451 shows an ex- 
tended 4a detection near the x-ray centroid (Donahue 1996), 
although the brightest cluster galaxy is offset from this po- 
sition by 23 arcsec, and the only radio source in the vicinity 
is offset to the other side by 20 arcsecs. Radio maps of the 
4 cluster centres with possible BCG detection, along with 
10 arcsec sub-mm error circles for reference, are shown in 
Fig. 4. 

Molecular gas, traced by CO emission, might also be 
present in strong cooling flows (Braine et al. 1995). For 
MS 1455, we also obtained an upper limit on the CO(4-3) 
emission from the BCG, which suggests a gas-to-dust ratio 
of ^ 150 if the SCUBA peak were actually a detection. This 
would be somewhat unusual for known properties of dusty 
BCGs (Edge et al. 1999), but a weak CO(l-O) detection for 
the MS 1455 BCG (A. Edge, private communication) lends 
support to a possible SCUBA detection. 

If the SCUBA detection in MS 0451 is associated with 
a cooling flow, then it is fairly unusual, and may indicate a 
massive, but highly obscured optical counterpart is actually 
being built by cooling flow fallout at the true centre of the 
cluster potential. 

We can place limits on the dust mass for the BCG 
for all clusters in our sample, by considering the 850 /xm 
flux density limit at the redshift of the BCG (see Edge et 
al. 1999), and adopting model dust parameters, for exam- 
ple Td =40K, with an emissivity index of 1.5. These limits 
along with the position of the BCG are presented in Ta- 
ble 3. Corresponding CO(l-O) detections and limits, as well 
as broader implications for cooling flows are presented in 
Edge et al. (in preparation). The limits are generally much 
less than the total mass deposition implied by the cluster 
cooling flows. But our observations are restricted to emis- 
sion around the BCG, and in any case the limits are model 
dependent. Any conclusions we draw must necessarily not 
be very strong. Our non-detections or relatively weak sub- 
mm emission from most BCGs in a sample of nine clusters, 
imply that substantial effort will be required to detect other 
central cluster galaxies. And that if some of the mass is being 
deposited as cool dust, this will be hard to detect in typical 
central cluster galaxies with current instrumentation. 
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Figure 4. VLA 4.9 GHz maps surrounding the briglitcst cluster galaxies in Zw3146, MS 1455, A 2261, and A 2219. The field size 
is 21 X 21 arcsecs in all cases. Radio contours are 3, 5, 10, 20, 50 and 100 times the RMS of the VLA map region. A 10 arcsec 
diameter error circle is overlaid to indicate the SCUBA beam position, for these marginal 2-3(t 850 /im detections. 









Table 3 










Limits 


on dust mass for central 


cooling flows. 




Cluster 


R.A." 


Dec." 


•=850 


Sj.n 


Dust mass^ 


Mo/yrd 




(.72000) 


(.12000) 


(mjy) 


(mJy) 


''^Mq 




CL 0016+16 


00 18 33.2 


+ 16 26 IS 


< 6.1 


< 66 


< 4.2 X 10** 


n/a 


MS 0451-03 


04 54 13.4 


-03 01 47 


19.1+4.2 


<354 


1.3 X uf 


n/a 


Abell520 


04 54 19.0 


+02 56 49 


<10.0 


<116 


< 2.5 X 10** 


<85 


Zwicky 3146 


10 23 39.6 


+04 11 10 


6.6+2.6 


< 48 


2.4 X 10** 


1358 


MS 1054-03 


10 56 06.2 


-03 37 27 


< 8.4 


<203 


< 8.5 X 10** 


n/a 


MS 1455+22 


14 57 15.1 


+22 20 35 


5.3+2.1 


< 42 


1.5 X 10** 


1227 


Abcll2163 


16 15 34.1 


-06 07 26 


<11.3 


< 86 


< 3.1 X 10** 


<90 


Abcll2219 


16 40 20.5 


+46 42 59 


3.6+1.7 


26 + 8 


9.0 X lO'' 


<485 


Abcll2261 


17 22 24.1 


+32 07 45 


9.5+3.8 


<205 


1.8 X 10** 


680 



'^ BCG position. 

^ Where no positive value > 2. Oct is detected, we quote Bayesian 95 per cent upper limits. 

== Assuming Td=40K and 13 = 1.5, with ^50 = //o/lOOkms^^Mpc^^ 

■^ Mass deposition values from Allen (2000). The three high- 2: clusters do not have the x-ray spatial resolution for accurate deprojection 

(Donahue et al. 1999) 



5 SOURCE COUNTS 

Source count models for rich cluster fields in the sub-mm 
have been studied in detail by Blain (1997) and Blain et 
al. (1999), with observational constraints presented in Small 
et al. (1997) and Ivison et al. (2000). For unmagnified blank 
fields, the sub-mm galaxy counts can be derived by simply 



dividing the number of detected sources by the surveyed 
area, although the effects of clustering may have to be con- 
sidered in small fields (Scott & White 1999). The situation 
is more complicated in the case of lensing due to a massive 
cluster, since the source plane is both distorted and magni- 
fied at a level dictated by the source redshift and position. 
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Thus in the image plane, various regions are observed to 
different depths. Although our clusters were not chosen to 
provide the largest lensing amplification, the models of Blain 
(1997, figure 4) show that our cluster redshifts are essentially 
optimal for maximum amplification, with the exception of 
MS 1054 which has slightly too high a redshift, but still in- 
curs a reasonably large lensing amplification. 

We model the amplification for each source and the ef- 
fective depth of the source plane area observed using an el- 
liptical potential model of each cluster (Korman et al. 1994, 
Newbury 1997). Parameters for the models are derived from 
both our own imaging data and from the best current infor- 
mation available in the literature. The lensing inversion is 
based on the prescription of Korman et al. (1994) and New- 
bury & Fahlman (1999), allowing a calculation of the lensing 
amplification and distortion as a function of position in the 
image plane. Table 4 presents the parameters for each clus- 
ter along with an estimate for the lensing amplification at 
each source position. Columns list the cluster name, redshift, 
cluster velocity dispersion, angular distance and position an- 
gle from North of each source relative to the x-ray centroid, 
and resulting lens amplification. The final column lists the 
reference for the cluster parameters. All distances assume a 
flat A-dominated cosmology with Q.m = 0.3. 

The essential parameters of the model are the redshift 
of the cluster and the 1-D (line-of-sight) velocity dispersion 
(cry), the latter being essential for setting the mass scale. 
These two quantities set the length scale in the image plane. 
Larger cluster redshifts reduce the effective magnification, 
while larger cluster velocity dispersions (cluster mass) in- 
crease the magnification. Although more detailed mass mod- 
els can in principle be constructed using the radial velocities 
of the individual cluster members (see Blain et al. 1999a), 
the errors in the source counts are of order 30 per cent, and 
our results would be very little affected by the use of such 
cluster mass models. 

There is still an ambiguity associated with the unknown 
background source redshift. However, the models can be 
used to constrain the range of possible magnifications. For 
sources at reasonably large redshift {z>l say), we find that 
the cumulative source counts depend fairly weakly on the ac- 
tual source redshift (less than 30 per cent scatter) in agree- 
ment with Blain et al. (1999a). If all sources are assumed 
to lie at z > 2, the scatter reduces further. We assume for 
simplicity that all sources lie at 2 = 2. Based on our CY 
radio analysis of the source redshifts, this is a reasonable 
assumption, although some sources may in fact lie at some- 
what lower redshifts, especially if their dust temperatures 
are colder than 50 K. 

The cluster ellipticity {E) is chosen based on an outer 
contour for the central gala^xy, if one exists. In the absence 
of this BCG contour, we assume circular symmetry. The el- 
lipticity is important only in the case of sources near the 
cluster centre such as in Zw3146. The core radius (i?coro) of 
the cluster is a model-dependent parameter - the core ra- 
dius of the x-ray gas being temperature dependent. However 
for clusters with lensed arcs in their image, it is possible to 
estimate the core radius directly. Again, the magnification 
is particularly sensitive to this number only near the cluster 
centre, as it sets the location of the critical lines with very 
large magnification. The ellipticity and core radius are quasi- 
independent parameters of the models. For modest values of 
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Figure 5. Amplification vs area in the source plane for the 8 
clusters in our sample with detected sources. 



the ellipticity, this parameter plays a very small role, espe- 
cially in view of all the other uncertainties involved. 

Lensing increases the brightness of sources and also 
magnifies the area - so we reach to fainter effective flux lim- 
its, but over a smaller effective solid angle. Fig. 5 shows the 
resulting magnification vs source plane area for each cluster. 
Our complete survey area in the source plane is in fact only 
~ 19 arcmin^, from a total surveyed SCUBA field area of 
~47 arcmin^. We estimate that the uncertainties in cluster 
parameters {ov, -Rcoro, E) to the lensing model contribute 
~ 25 per cent error, a little larger than the errors in cal- 
ibrating SCUBA maps, but less than the Poisson error in 
the source counts themselves. 



6 RESULTS 

With the lensing factored in, our effective noise level in 
the central regions of the cluster fields decreases sharply, 
as displayed in Fig. 6. From these relations of sensitivity 
as a function of cluster area, we can calculate the effective 
area surveyed to a given depth. Fig. 7 shows the integral 
source counts from our survey including all sources detected 
at > 3a (solid circles, binned). We plot the counts from the 
Small, Ivison, Blain (1997) lensing survey presented in Blain 
et al. (1999c, crosses), as well as a fit to the combined blank 
field sub- mm counts (Hughes et al. 1998, Barger et al. 1999b, 
Bales et al. 1999) constrained at faint fluxes to not overpro- 
duce the far-IR background (dotted line). When we restrict 
ourselves to the > 4cr sub-set which contains no obvious 
cluster members there is no dramatic change in the counts 
(open circles in Fig. 7). We fit these 4cr counts with a sim- 
ple power law A''(>S') = A''o(S/So)~^, with So = 1 mjy, and 
find 7~ 1.73, and A^o — 11200. These are close to the blank 
field values of 7 = 2.2, and No = 13600. A direct comparison 
with other published sub-mm source counts thus reveals a 
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Tabic 4 
Lcnsing parameters for clusters 



Cluster 


z 


(Tv{krn s 






A9sourco (PA)" 






Lens 


amplification 


Cluster reference 


CI 0016 + 16 


0.541 


1243 






44(138), 62(171) 








3.0, 2.6 


Ellingson et al. 1998 


MS 0451-03 


0.550 


1353 






65(130), 25(46) 








2.6, 4.5 


Donahue (1996) 


Abell520 


0.203 


927 






61(104), 74(138) 








3.4, 3.1 


Gioia & Luppino (1994) 


Zwicky3146 


0.291 


1310 


31 


117) 


62(179), 69(-33), 


78( 


136) 


6.8 


3.7, 3.5, 3.2 


Ebeling et al. (1996) 


MS 1054-03 


0.833 


1170 






63(70) 








2.6 


Donahue et al. (1998) 


MS 1455+22 


0.259 


1168 






58(162.5) 








2.8 


Stocke et al. (1991) 


Abell2163 


0.201 


1231 






n/a 








n/a 


Gioia & Luppino (1994) 


Abell2219 


0.228 


1000 


66(17), 


45(157), 15(178), 


41(-93) 


3.6 


3.2, 7.5, 3.0 


Small ct al. (1996) 


Abell2261 


0.224 


1102 






72(-124) 








3.3 


Ebeling ct al. (1996) 



Distance and position angle of the sources (listed in Table 2) from the x-ray centre. 
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Figure 6. Area vs sensitivity for the 8 clusters in our sample 
with detected sources. Line types are the same as in Fig. 5. 



good agreement between results, including gravitationally 
lensed and blank field counts. SCUBA preferentially selects 
galaxies with z ^ 1, due to the spectral shape at these wave- 
lengths. Therefore we expect few detections in the clusters 
themselves, however, this may be the source of some of the 
small discrepancy between the various surveys. 

A positive magnification from the cluster lensing means 
that we are seeing fainter galaxies, and hence an increase 
in the surface density of galaxies. However, this is offset by 
the loss in area due to field distortion. If the source count 
power-law index were too shallow, this would result in a 
less efficient survey compared to a blank field survey with 
equivalent integration time. However, our results show a net 
positive bias, so that we are detecting more sources per unit 
time and at unlensed flux densities below what would be 
possible in blank field searches, in line with what was found 
earlier by Smail et al. (1997). 

Using the amplification due to the foreground cluster, 
lensed surveys also suffer less from confusion noise (e.g. Blain 
et al. 1998). The large beamsize of SCUBA implies a con- 
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Figure 7. The integral source counts (number of objects brighter 
than flux density 5 as a function of S. We plot our total counts 
of > 3(T detections (solid circles, binned), and our sub-sample ex- 
cluding possible cluster members and detections less than 4cr sig- 
nificance (open circles, individual sources). The corrected counts 
from Blain et al. (1999c) are plotted with crosses. The dotted line 
depicts a fit to the combined blank field sub-mm counts (Hughes 
et al. 1998, Barger et al. 1999b, Bales et al. 1999) constrained at 
faint fluxes to not overproduce the far-IR background. The error 
bars shown for our complete sample points are ±1ct (68 per cent 
Bayesian confldencc region) errors based on Poisson statistics. 



fusion limit at ~ 1.5 mJy, while we have detected 4 possible 
sources fainter than 2mJy. 



7 DISCUSSION 

Estimates of the redshift distribution for the SCUBA- 
selected sources remain in dispute (Hughes et al. 1998; 
Barger et al. 1999a; Lilly et al. 1999a). Hughes et al. (1998) 
concluded that the bulk of the population is at 2 --^ 2- 
4, based on photometric redshift limits for the probable 
counterparts of five sub-mm sources in the Hubble Deep 
Field (c.f. Richards 1999; Downes et al. 1999). Barger et 
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al. (1999a) undertook a spectroscopic survey of the Smail et 
al. (1999) sub- mm sample and concluded that the median 
redshift was ~ 1.5-2, with the bulk of the population hav- 
ing 2~l-3. Lilly et al. (1999a) used archival spectroscopy 
and broad-band photometry of sub-mm sources from the 
Eales et al. (1999) survey to conclude that the population 
spans 2 = 0.1-3, with a third at z<l. Given the relatively 
small number statistics and the very real possibility of some 
misidentifications, it is unclear whether there is any genuine 
disagreement. 

However, use of the Carilli & Yun (CY - 2000) rela- 
tion of radio/far-IR flux to predict the redshifts for a large 
sample of SCUBA sources (Smail et al. 2000), has suggested 
that the population may have a median value lying between 
2 = 2.5 and 3. This is consistent with our CY analysis using 
VLA maps, which result in a median value lying at 2 <; 2.0, 
although dependent on the dust temperature and emissivity 
adopted. Fig. 8 shows an estimate of the cumulative red- 
shift distribution for our entire sample, derived from the 
CY relation using radio detections and lower limits. Also de- 
picted are the cumulative redshift distributions for models 
formulated to fit the far-IR background, and various infrared 
counts, by Guiderdoni et al. (1998), and Blain et al. (1999b). 

The archival radio maps for many of our cluster fields 
are not deep enough to put strong constraints on the red- 
shift distribution of sub-mm sources, with few actual ra- 
dio detections. Given that the derived redshifts are gen- 
erally lower limits, it is difficult to presently discriminate 
between predictions for the redshift distribution from var- 
ious models (e.g. Blain et al. (1999b) model Gaussian and 
modified Gaussian models (long and short dashes) or the 
Guiderdoni et al. (1998) model (dotted line)). However, the 
radio/far-IR correlation indicates that the bulk of the dust 
emission in the universe did not occur at redshifts much 
below 2 = 2. The true redshift distribution remains to be 
determined. Deep optical and near-IR imaging should also 
uncover these objects' optical properties, and resolve cur- 
rently unanswered questions such as whether the galaxies 
would be detectable through their Lyman break (Chapman 
et al. 2000), or whether they have extremely red colours 
(Smail et al. 1999). 

We can calculate lower limits to the far-IR background 
radiation intensities from measured flux densities of resolved 
4cr sources in our cluster fields. The contribution is the sum 
'E{S/A), where A is the area which has a sensitivity neces- 
sary to see a source of flux S, where lensing has been factored 
into the sum. We find 9.1 Jydeg"^ averaged over all 9 fields. 
This corresponds to 20 per cent of the total FIRB at 850 /im 
as measured by Fixsen et al. (1998): 44Jydeg^^. The ac- 
curacy of these estimates however depends on the lensing 
models and other details of the counts, together with the 
remaining uncertainty in the FIRB measurement as these 
wavelengths. If we take the simple power law model inte- 
grated up from 1 mjy, then the background is A^o7/(7 — 1). 
Using our slope from Fig. 7 accounts for more than half of 
the background down to 1 mJy. 
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Figure 8. Cumulative redshift distribution estimated for our en- 
tire sample, derived from the CY relation using VLA data (solid 
line). Also shown are the cumulative redshift distributions for 
models formulated to fit the far-IR background, and various in- 
frared counts, by Guiderdoni et al. (1998) — dotted line, Blain et 
al. (1999b) Gaussian model of luminosity evolution - long-dashed 
line, and Blain et al. (1999b) modified Gaussian model fit to the 
Barger et al. (1999a) redshift distribution - short-dashed line. 



suggests that most of these are being lensed by the cluster 
mass, rather than being in the clusters themselves. 

(ii) We see no strong evidence for cool dust emission 
around BCGs, which might be related to cooling flows. 

(iii) Our sub-mm source counts are in reasonable agree- 
ment with existing sub-mm surveys, and extend the list of 
such sources available for follow up at other wavelengths. 

(iv) The redshift distribution has been loosely con- 
strained using the radio- far-IR correlation, and favours a 
relatively high redshift (2 > 2) for the bulk of dust emission 
in the Universe. 
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8 CONCLUSIONS 

(i) We have detected 17 new sub-mm sources in the fields 
of 9 rich galaxy clusters. Comparison with VLA radio images 
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